import control as ctrl   #pip install control
import matplotlib.pyplot as plt
import math
import numpy as np
sysTfList=[]
sysTfList.append(ctrl.tf([1],[1,0.4,0]))
mag,phase,omega=ctrl.bode(sysTfList,omega=np.logspace(start=-2,stop=1,num=200),dB=True,deg=True)

plt.show()